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The first-order correction of the perturbative solution of the coupled equations of the quadratic 
gravity and nonlinear electrodynamics is constructed, with the zeroth-order solution coinciding with 
the ones given by Ayon-Beato and Garcia and by Bronnikov. It is shown that a simple generalization 
of the Bronnikov's electromagnetic Lagrangian leads to the solution expressible in terms of the 
polylogarithm functions. The solution is parametrized by two integration constants and depends 
on two free parameters. By the boundary conditions the integration constants are related to the 
charge and total mass of the system as seen by a distant observer, whereas the free parameters are 
adjusted to make the resultant line element regular at the center. It is argued that various curvature 
invariants are also regular there that strongly suggests the regularity of the spacetime. Despite the 
f"^ ■ complexity of the problem the obtained solution can be studied analytically. The location of the 

£SJ ' event horizon of the black hole, its asymptotics and temperature are calculated. Special emphasis 

is put on the extremal configuration. 
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t-H , I. INTRODUCTION 

\Q . 

One of the most important and intriguing questions of the black hole physics is the problem of singularities that 
reside in their internal region, hidden to an external observer by the event horizon. In the vast majority of papers 
, singularities are treated as symptoms of illness of the theory rather than its health (see however Ref. [jj for a different 
point of view), and, consequently, a great deal of efforts were directed to constructing singularity-free models. However, 
a subtle point is that the Einstein field equations loose their predicative power and cannot be trusted when the 
curvature of the manifold approaches the Planck regime. Indeed, according to our present understanding a proper 
description of the gravitational phenomena should be given by the quantum gravity, being perhaps a part of a more 

• • ■ fundamental theory. And although at the present stage we have no clear idea how this theory looks like, we expect 
that the action functional describing its low-energy approximation should consist of higher order terms constructed 
from the curvature and its covariant derivatives to some required order. It means that in the full theory the analogs 
of the nonsingular solution of the Einstein gravity may loose their nonsingular character as well as the singular ones 

5^ i their singularities. 

Among various modifications of the general relativity the prominent role is played by the quadratic gravity 0, 
H, 0, IE IE 111 Motivations for introducing into the action functional terms which are quadratic in curvature are 
numerous. For example, when invented, the equations of quadratic gravity have been treated as an exact formulation 
of the theory of gravitation. For historic informations and important references the reader is referred to Ref. [9(. It 
may be considered, quite naturally, as truncation of series expansion of the action of the more general theory. Such 
terms appear generically in one- loop calculations of the quantum field theory in curved background |lOj| . Moreover, 
from the point of view of the semi-classical gravity it might be treated (in certain circumstances) as some sort of a 
poor man's stress- energy tensor, allowing in a relatively simple way to mimic, especially when the application of the 
full stress- energy tensor would produce extremely complicated or even intractable results, the fairly more complex 
source term of the field equations. 

Analyses of the spherically-symmetric and static solutions to the higher derivative theory has been carried out in 
ITTl u3. IT3I. E3 . Il5 | . For example, in Ref Q it has been shown that the weak- field limit of the quadratic gravity 
involves, beside the ordinary Newtonian term, also the terms with the Yukawa-like potential. Series solutions near 
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the r = have been investigated in [3|, LLJ] ■ Such solutions are limited to the closest neighborhood of the center and 
should be matched to the appropriate solutions valid for large r that require numerical integration. However, an 
important lesson that follows from this calculations is that the regular solutions of the equations of the higher order 
theories are quite common 

As compared with the General Relativity, the Lagrangian of the quadratic gravity in four dimensions requires two 
additional terms aR a bR ab + flR 2 , where a and j3 are the coupling constants. Since the coupling constants are expected 
to be very small, one can easily device a perturbative approach to the problem treating the classical solution of the 
Einstein field equations as the zeroth-order of the approximation. Successive perturbations are therefore solutions 
of the differential equations of ascending complexity. It should be noted that although the method is clear the 
calculations beyond the first- order may be intractable. 

According to the well-known theorem, if the Lagrangian C (F) has Maxwell asymptotics for weak fields, i.e., C (F) ~ 
F and £f = dC/dF — ► 1 as F — > 0, then any static and spherically-symmetric solutions to the system of coupled 
equations of nonlinear electrodynamics and general relativity characterized by electric charge, Q e , cannot have a 
regular center [Ifil IrH ITsj. This no-go theorem does not forbid existence of the solutions with magnetic charge, Q, 
as well as some hybrid solutions in which the electric field does not extend to the central region. On the other hand, 
however, it has been argued in Ref. [19j that the Maxwell asymptotics at great distances are essentially different from 
these at the center, and, consequently, the condition Cf — ► 1 as _F — > is too restrictive. 

The issue of the regular black holes in general relativity has a long and interesting history. For example, one of the 
method that can be used in construction of such configurations is simply re placing the singular black hole interior by 
a regular core. This idea appeared almost forty years ago, in mid sixties [2(j, |2l], |22] and is actively investigated today. 
In models considered in Refs. [2^, [24| part of the region inside the event horizon is joined through a thin boundary 
layer to de Sitter geometry. Such a geometric surgery certainly does not exhaust all interesting possibilities: of equal 
importance are the regular geometries with suitably chosen profile functions, or, better, exact solutions constructed 
for specific, physically reasonable sources [2^, |2(| |23, IH, |2!j |3(j • One of the most intriguing regular solutions to the 
coupled equations of the nonlinear electrodynamics and gravity have been constructed by Ayon-Beato and Garcia j^H 
and by Bronnikov ^tJ- We shall refer to solutions of this type as ABGB geometry. The former describes a regular 
static and spherically symmetric configuration with an electric charge whereas the latter describes a similar geometry 
characterized by the total mass M and Q. For certain values of the parameters both solutions describe a black hole. 
It should be noted that the electric solution does not contradict the non existence theorem as the formulation of 
nonlinear electrodynamics employed by Ayon-Beato and Garcia (V framework in the nomenclature of Refs. |171 Il8|') 
is not the one to which one refers in the proviso of the no-go theorem. Indeed, the solution of Ayon-Beato and Garcia 
has been constructed in a formulation of the nonlinear electrodynamics obtained from the original one (T framework) 
by means of a Legendre transformation (see Ref. 01 f° r details) . An attractive feature of this solutions that certainly 
sim plifi es calculations is possibility to express the location of the horizons in terms of the known special functions 

mm. 

The natural question that arises in connection with the foregoing discussion is whether or not it is possible to con- 
struct an analog of the solution of the ABGB-type in the quadratic gravity which shares with its classical counterpart 
regularity at r = 0. And although the full, detailed answer is beyond our capabilities, it is possible to provide an 
affirmative answer to the restricted problem. Indeed, since the complexity of the coupled equations of the quadratic 
gravity and nonlinear electrodynamics, even in the simplest case of spherically-symmetric and static geometries, 
hinders construction of the exact solution, one has to refer to the analytical approximations or numerical methods. 

In this paper we shall employ perturbative methods to construct the first-order solution to the equations of the 
quadratic gravity with the source term being generalization of the stress-energy tensor of the Bronnikov type. The 
zeroth-order solution coincides, as expected, with the ABGB line element whereas the first-order correction can be 
elegantly expressed in terms of the polylogarithm functions. An interesting feature of this very solution is its regularity 
for r — 0. The Kretschmann scalar and other curvature invariants are also regular at the center that suggests regularity 
of the underlying geometry. It should be emphasized that the demonstration of the regularity of the full solution 
would require either profound understanding of the perturbation series to any required order or construction of the 
full, physically acceptable nonperturbative solution. On the other hand, the perturbative approach is expected to 
yield reasonable results especially for higher derivative dynamical equations. In fact it may be the only method to deal 
with them. Indeed, since the quadratic gravity involves fourth-order derivatives of the metric their nonperturbative 
solutions may appear to be spurious and one has to invent a method for systematic selecting physical ones. It seems 
that the acceptable solutions, when expanded in p owers of the small parameter, should reduce to those obtained 
within the framework of perturbative approach |34l I3M l3( | . 

The paper is organized as follows. In Sec. [H] we introduce basic equations and briefly sketch employed method. We 
choose the line element in the form propounded by Visser [3?|. which has proved to be a very useful representation 
considerably simplifying calculations. In Sec. II I II we introduce the Lagrangian of the nonlinear electrodynamics, 
construct solutions to the first-order equations and establish regularity of the thus obtained line element. Subsequently 
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we study the limit of the vanishing Q and analyse the regularity of various curvature invariants. Corrections to the 
location of the inner and outer horizons and to temperature are given in Sec. II V Al The position of the horizons of 
the ABGB spacetime is given in terms of the real branches of the Lambert functions. The extremal configuration 
is studied in Sec. IIV Bl Specifically, we calculate modifications of the location of the degenerate horizon caused by 
quadratic gravity and analyse relations between Q and the total mass as seen by a distant observer. Sec. IV1 contains 
final remarks. In Appendix we have collected useful formulas and presented a brief description of the method of 
integration of the field equations in terms of the poly logarithm functions. Throughout the paper the geometric 
system of units has been adopted and the the sign conventions are taken to be that of MTW |33 ■ 



II. THE EQUATIONS 

In absence of the cosmological term, the coupled system of the nonlinear electrodynamics and the quadratic gravity 
is described by the action 

S = — — ~p;S g + SVra, (1) 



16ttG 



where 



and 



S g = J (R + aR 2 + l3R ab R ab ) V=9~d 4 x, (2) 



S m = --^- I £(F)V=gd 4 x. (3) 



16tt. 

Here C (F) is some functional of F = F a i,F ab (its exact form will be given later) and all symbols have their usual 
meaning. The third possible term constructed form the Kretschmann scalar, r yR a bcdR abcd , may by removed from the 
Lagrangian with the help of the Gauss-Bonnet invariant. Of numerical parameters a and (3 we assume, as usual, 
that they are small and of comparable order, otherwise they would lead to the observational consequences within 
our solar system. Their ultimate values should be determined form observations of light deflection, binary pulsars 
and cosmological data 

SEE! 

Moreover, following Ref. 01, we shall restrict ourselves to spacetimes of small 

curvatures, for which the conditions 

\aR\ « 1, \/3R ab \ « 1 (4) 

hold. Although demanding that the mass scales associated with the linearized equations are real may place additional 
constraints |42j, |43J, |4J] on a and /3, we shall treat them as small but arbitrary. 
The tensor F a b satisfies equations 

v „(^m^)=o, (5 ) 

V a *F Qb = 0, (6) 
and the asterix denotes the Hodge duality. The stress-energy tensor defined as 

T ab = * 6 S m (7) 

V~9 °9ab 



is given therefore by 



Differentiating functionally the action S with respect to the metric tensor one has 

L ab _ Qab _ a jab _ pjab = g^a^ (g) 
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where 



and 



jab = m ab _ 2RR ab + }_ g ab [ R 2 _ ^ 



jab = R ,ab_ aR ab _ 2 R cd R cMa + \g ah (R cd R cd ~ OR) . (11) 

Let us consider the spherically symmetric and static configuration described by the line element of the form 

ds 2 = -e 2 ^f(r)dt 2 + ^-+ rW, (12) 



where 



2M(r) , x 

/(r) = l-— (13) 



The spherical symmetry places restrictions on the components of F a b tensor and its only nonvanishing components 
compatible with the assumed symmetry are Fqi and 1*23- Simple calculations yield 

F 23 = Qsm0 (14) 

and 

re — F w =Q e , (15) 

where Q and Q e are the integration constants interpreted as the magnetic and electric charge, respectively. In the 
latter we shall assume that the electric charge vanishes, and, consequently, F is given by 

F=%. (16) 
The stress-energy tensor (JJJ calculated for this configuration is 

Tl=Tl = -±-L{F) (17) 

and 

To simplify calculations and to keep control of the order of terms in complicated series expansions we shall introduce 
a dimensionless parameter e substituting a — > ea and /3 — > e/3. We shall put e = 1 at the final stage of calculations. 
Of functions M (r) and ip (r) we assume that they can be expanded as 

M (r) = M (r) + eM x if) + O (e 2 ) (19) 

and 

V> (r) = # a (r) + O (e 2 ) . (20) 

Consider the left hand side of Eq. © calculated for the line element (|12|l first. Making use of the above expansions 
and collecting the terms with the like power one obtains 

L\ = -^{M' +sM[-eSl), (21) 
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where 



S$=f3 



2M; 8M m; 2M; 2 2M; 5M M£ Mq Mq 



+M! + M f - _ Mg M (3) + r Af < 4) - 2 Af M< 4) j 

_ a / 24A/ A^ _ 8MJ _ 4Af^ 2 + 8A^ _ 18M Afff _ M „ 2 



« _ 4 M (3) + 6MoMf + 2 M , M ( 3) _ 2 r M (4) + 4 Mo M ( 4 



(22) 



and Mq, Mq and Mq for i > 3 denote first, second and i— th derivatives with respect to the radial coordinate. On 
the other hand, a simple combination of the components of L b a tensor 

L r r -L\=Q (23) 

can be easily integrated to yield 

Mr) = (2a + /3)M< 3) - 4r(3a + (3)M' Q + C u (24) 

where C\ is the integration constant. It should be noted that contrary to the case of coupled system of the Maxwell 
equations and quadratic gravity considered in Refs. [TlL 1 1 2l f]~^ . I X f5l | now we have explicit dependence on the parameter 
a. This together with the nonlinear character of the source term (that will be specified below) results in substantial 
complications of the first-order equations. 

To develop the model further we shall express solutions of the system of differential equations consisting of the time 
component of Eqs. © and Eq. l|23|l in terms of the total mass M. as seen by a distant observer 

lim A/(r) =M, (25) 

r — >oo 

whereas for the function if> (r) we shall adopt the natural condition 

lim V (r) = 0. (26) 



III. SOLUTIONS 



Further considerations require specification of the Lagrangian C (F) . Following Ayon-Beato, Garcia and Bronnikov 
let us chose it in the form 



C(F) =F 




1-tanhM,^ 



(27) 



where 



26 ' 



(28) 



and the free parameter b will be adjusted to guarantee regularity at the center. 

Before proceeding further we shall briefly discuss the question of regularity, postponing presentation of the technical 
details for a while. First, observe that the zeroth-order solution coincides with a general ABGB line clement that 
depends on a free parameter b and an integration constant C2 ■ On the other hand, the first-order solution written in 
the suitable representation approaches, as we shall see, a constant value, say fj,, as r — > 00, and, consequently, by the 
boundary conditions (|25|l one has C2 = M.—e[i. The thus obtained line element is generally singular at the center, and 
the only way to make it regular consists in appropriate choice of the free parameter. The regularity is understood here 
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as the regularity of the line element rather than regularity of the spacetime itself as the latter requires the curvature 
invariants to be regular. We shall return to this issue later. 

Now we present the calculations in a more systematic form. In order to guarantee sufficient generality of our 
considerations we shall take the parameter b in the form 

b = b 1 +eb 2 , (29) 

with b 2 7^ 0. Since we have assumed the expansions of M (r) and ip (r) in the form given by Eqs. (|19|l and 12Ufl . 
respectively, we shall rewrite the boundary conditions as 

lim M (r) =M, (30) 

r — >oc 

with vanishing Mi(oo) and ip^oo). 

Inserting Eq. H28|) into l)27|) and makig use of Eq. one has 

/W-2£(i-t-tf£). (3D 

Subsequently expanding the right hand side of Eqs. @ with respect to s and retaining the linear terms only yields 



8nT* = 8ttT; = -^L. 1 - tanh 2 ^— - e^V cosh^ 2 tanh ^— . (32) 
11 2bxr) b\r° \ 2b x r 2hr } 



The zeroth-order equation 
can be easily integrated 

M (r) =C 2 -6 1 tanh^-, (33) 

where C 2 is the integration constant. Making use of the condition (|30ll gives C 2 = M. On the other hand, demanding 
of the regularity of the line element as r — > yields b\ — M , and, consequently, Mq (r) reads 

M (r)=M(l-tanh^y . (34) 

The zeroth-order solution has an interesting property, which, as we shall see, is crucial in our subsequent analysis: 
Mq(t) as well as its derivatives vanish in the limit r — > 0. Moreover, it should be emphasized that although the 
integration constant and the free parameter in the thus constructed solution are equal to the total mass of the system 
as seen by a distant observer their status is different: the interpretation of the former is an inevitable consequence of 
the boundary conditions whereas the latter should be postulated. 

For small values of |Q|/.M as well as at great distances the ABGB line element resembles that of Reissner-Nordstrom. 
It can be easily seen by expanding the function Mo(r) 

2M Q 2 Q 6 , 

9tt = "I + 2" + 4 + -' ( 35 ) 

r 12.M r 4 



AM f-Q 2 , 

/ ~ 1 exp —— . (36) 



whereas for r — > one has 

L 

r \ Mr 

Noticeable differences appear for the configurations near the extremality limit and in the internal region in the vicinity 
of the center. Indeed, for the ABGB solution g tt tends to — 1 as r — > whereas the (00) component of the metric 
tensor of the Reissner-Nordstrom solution diverges in that limit as — r~ 2 . 

From igU one concludes that the zeroth-order solution, Mq (r), is sufficient to determine the function ip 1 (r). 
Indeed, substituting l|34|) into Eq. I|24|) and making use of the condition (|26|l one obtains 



1>i ( r ) 



PQ 2 , {2a + 0)Q\ ^ 2 Q 2 3(2a + p)Q\ u Q 2 



tanh — — — -= tanh ■ 



r 4 2A^ 2 r 6 2Mr Mr 5 2Mr 

(2a + f3)Q\< g : 



cosh 2 



2Mr 



cosh- 4 -f— . (37) 



AM 2 r e 2Mr 
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It could be easily demonstrated that that ipi(r) vanishes at r = and inspection of Eq. I|24|) reveals similar behaviour 
of its derivatives. 

The solution of the first-order equation 



A(M{-^) = |^fcosh- 2 -- 
r l b\r n \ 2b\r 



Q" , Q 2 

tann 



2b ir 



is more complicated and when combined with the appropriate boundary conditions it could be written as 



Mr (r) 



Sl(r)dr + 



b 2 Q 4 
2M 2 



~% ( cosh > i « 
ri \ 2Mr 



tanh - 



2Mr 



dr. 



(38) 



(39) 



Let us consider the first integral in the right hand side of l|39|l first. The result can be expressed in terms of the 
hyperbolic functions and polylogarithms Li, (s). Indeed, utilizing formulas collected in Appendix one can construct 
the solution which has the general structure 



J oo •' Vl • r> .-_n 



i=0 j=0 



--^£(^ a) +/^ Lii (" exp 

4 — 1 ^ 



.Mr 



2 Air 



2 Air 



(40) 



Here // is given by 



TT 2 M 5 



31 



315 



31 



— — 7T + /3 - + — 71^ - — 7T 



45 



630 



7T 2 >1 5 
Q 6 



(41) 



and since the terms in the square brackets will appear frequently in our subsequent analyses we have singled them out 
and denoted by a. The functions / and h are simple polynomials of r _1 ; their actual form will not be displayed here 
as we shall readily rewrite them in a slightly modified form. It should be noted, however, that a careful term-by-term 
analysis of Eq. 1)40(1 reveals that M{ (r) approaches —fi as r — > 0. Since for b — b± the functions M\(r) and 
coincide it is impossible to obtain a regular solution at the center, and the remedy is to retain b 2 in the calculations. 
Converting all hyperbolic functions appearing in the solution to the exponents, introducing a new function £ defined 

as 



f (r) = exp - 



Mr 



(42) 



and finally making use of the elementary properties of the function Lio(— £) (see Eq. (|A.11|) of the Appendix) one has 



M 



i=0 



(43) 



The term fi has been singled out for convenience and and fj are given respectively by 



(a) 



48M 3 4MQ 2 (f - 1) 



16M' 



4Q 6 (1 - 4£ + T 



Q 2 r 2 r 4 (l + £) 2 r 3 (l + C) 7W 2 r 5 (l + 3 
4Q 4 [X (1 + 86£ - 89£ 2 ) + 25r (£ 2 - l)] 2Q 6 £ (75 - 425£ + 125£ 2 + f) 



5Mr 5 (l+0 



15Mr 6 (l + 



(44) 



M = 4Q 4 _ 96M 4 
Jl 5r 5 Q 4 r : 



{a) 4MQ 2 _ 96M 5 

h r 4 g 6 ' 



(45) 



6r 3 A4 3 (a) 
Q 6 /3 



2r 2 A1 2 (a) 



rM 



Q 2 h 



(a) 



f(«) 



96A4 £ 



(46) 
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and 



Jo 



(0) 



16M 3 _ 16M 2 _ 2Q 2 [Ur (1 + g) + M (3 - 36g + g 2 )] 
' Q 2 r 2 " 3r3(l + C) 3r 4 (l+0 2 
2Q 6 (l - 4£ + C 2 ) <3 6 £ (75 - 425£ + 125£ 2 + <f) 



7W 2 r 5 (l + 157Wr 6 (l + 

2Q 4 [M (1 + 96£ - 109£ 2 ) + 30r (<f - l)] 



(47) 



w 2Q 4 _ SAf_ _ 32A4 4 (/3) 



32.M 5 2A4Q 2 8X 3 



(48) 



m = 87W 2 _ 16M 4 



f(J3) _ _ 
■14 — 



16M 5 247W E 



Q6 g2 r 2 



(49) 



— f W = f0») = 1^ 
Q 2 1/5 /e Q6 



/f ) =/r = ™. (50) 



The second quadrature in the right hand side of Eq. I|39|l is elementary and gives 

— 2 .^o / -7 I cosh -2 J^— tanh \ dr = cosh -2 J^- b 2 tanh -rr— (51) 

1 2A4 2 J x r 3 V 2Mr 2Mr J 2Mr 2Mr 2Mr v ' 

or equivalently 

1 - £ 2Q 2 



£ Mr(l + 



Li (-0- (52) 



Now, in order to obtain a regular solution at the center one has to adjust suitably the parameter b 2 . To accomplish 
this let us observe that (0) = —b 2 . Further, taking the results of the discussion presented below Eq. (|41|l into 
account one finds 

b 2 = - M . (53) 

It could be easily seen that with such a choice of b 2 the function Mi (r) vanishes at the center and as r — > 00 as 
required. 

It should be noted that the representation of the thus obtained line element is by no means unique. For example 
one can express the result in terms of Li n (— 1/£) rather than Li n (— £) employing identity |45| 

[n/2] n _ 2r 

Li n (-£) + (-l)"Li n (-l/£) = ~-ln"g + 2 V D *, Li 2r (-l), (54) 

n\ ' (n — 2r 

where [n/2] is the greatest integer contained in n/2 and the constants Li 2r (— 1) are related to the Bernoulli numbers, 
B 2r , according to the formula 

Li2r(-1) - %^ B ^ 2r - (55) 

Regardless of the chosen representation, after imposing boundary conditions, both solutions are, of course, equivalent. 

As the ABGB geometry reduces to the Schwarzschild solution in the limit Q = 0, it is the charge, no matter how 
small, that secures regularity. The vanishing of the charge leads therefore to dramatic changes in the geometry of the 
black hole interior. On the other hand, the terms calculated to 0(Q 2 ) for the ABGB spacetime coincide with the 
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Reissner-Nordstrom solution. Let us consider a series expansion of the function m} 1 -* (r) for small q = |Q|/.M. After 
some algebra one has 



Mi = 



- A„ 2 

x 3 x i )M q 



-a 



<l 



-0 

7 ' ) Mx 7 



6/3 4 
~5M q 



+ ... - 6 2 



11 

2^ 



12a; 3 



60a; 5 



9 4 
4a; 5 a; 6 ' 



6720a; 7 / 



(56) 



where x = r/M. It could be easily checked that for the regular line element the leading term of the expansion is 
proportional to 62Q 6 . Since the coefficient 62 is proportional to Q~ e the constructed line element does not approach 
the Schwarzschild solution in the limit Q — > 0. Indeed, simple calculations yield 



ds 2 = - 1 - 



2M 2M 2 



k)dt 2 + 



i-^'Vi-^Y 



2M 2 



r J 



dr 2 +r 2 dn 2 , 



(57) 



where k is given by 



k = —a w 5.781a - 0.4865. 
12 H 



(58) 



Consequently, one has either Schwarzschild asymptotic of a singular line element with 62 = (fc = 0) or a regular 
line element with non-Schwarzschild Q = limit for 6 2 = — /1. It is possible, of course, to accept other values of the 
parameter 62 but it seems that they are of lesser importance. 

The approximate location of the event horizon of the line element (|57|) lies near its Schwarzschild value and is 
approximately given by 



r+ m 2M 1 + 



8M< 



(59) 



As is well-known, the Hawking temperature, Th, can be easily calculated from the metric tensor without referring to 
the field equations. The standard by now method of obtaining Tjj relies on the Wick rotation. The Euclidean line 
element has no conical singularity provided the time coordinate is periodic with a period P given by 



P = 4-7T lim {g t t9rr) 1/2 [ -j-gtt 

r^r + \ dr 



(60) 



Its reciprocal is identified with the black hole temperature, which, in the case in hand, reads 

1 / k 



T H = 



8nM 



1 



4M< 



(61) 



Note that for a and (3 satisfying 



1680 + 98tt 2 - 31tt 4 
2(31tt 4 - 2520) 



(62) 



both (|57|l and l|61|l reduce to their Schwarzschild counterparts. 

Finally, we shall investigate the important issue of regularity. However, before we proceed further a few words of 
comments are in order. It should be emphasized that we have constructed a linearized solution to the coupled system 
of equations of quadratic gravity and nonlinear electrodynamics only. Consequently, when calculating the curvature 
invariants we are to restrict ourselves to 0(e) terms. As it has been stressed earlier, the regularity of the line element 
does not necessarily entail the regularity of the underlying geometry: the latter requires various curvature invariants 
to be regular in the interesting region. First, let us concentrate on the Kretschmann scalar, K. It could be easily 
shown that for the constructed line element it consists of terms involving products of the functions Mo(r), Mi(r), 
if>i(r) and their derivatives. Indeed, after some algebra one obtains K = Kq + eAK, where Kq is the Kretschmann 
scalar of the ABGB spacetime 



Kn = 



48M n 2 



16M(5'Mo 



32M^ 



MM^Mq 



ami; 2 



(63) 
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and A if is given by 



fl28M M(y 16M 80Af 2 16M^ 24M M V ) / 24M[ ) M y 48M', 2 \ , 
AK = I ^ + — ~5 ^ ^ + ~2 — J Vi 

8M(! 16M M ' 32M 2 16M^ 32M M( ) 16M 



IGMq 96M 64M^\ /16MJ 16M 8M£'\ 

The regularity of K$ has been demonstrated in Ref |3l|; in order to establish regularity of AK one has to establish 
regularity of its building blocks. To analyse behaviour of the k— th derivative of M\ we shall employ the first order 
equation. From Eq. (|22|) one clearly sees that all derivatives of S\ vanish as r — > 0. It is simply because k— th derivative 
of Mq has the asymptotic form 



2k 

Ci 



)2 



-jpr~ E ^(-m;)' (65) 

i=k+l v 7 

as r — > 0, where the coefficients c, depend on M. and Q. It follows then that fc-th derivative behaves as 
r~ 2k exp(— Q 2 /Air) near the center. One encounters a similar behaviour in the second term in the right hand 
side of Eq. (|39|l . Pulling this all together one concludes that the Kretschmann scalar is indeed regular at the center. 
Moreover, one can easily establish the regularity of R 2 and R a i,R ab . Similar arguments can be used in demonstration 
that higher invariants constructed from the curvature are also regular. 



IV. GEOMETRY 



Generally speaking a black hole belongs to one of the two distinct classes: it may be either extremal (when horizon 
has at least twofold degeneracy) or nonextremal. In the former case the family of ultraextremal black holes, i. e. 
configurations with triple (or even higher) degeneracy can be singled out |4y, l47j . Since we have restricted ourselves 
to the case of vanishing cosmological constant, we expect that the ABGB black hole possesses at most two horizons. 
(Generalization of the ABGB solution to A 7^ case is straightforward.) Indeed, depending on q = \Q\/A4 the ABGB 
spacetime has two, one, or has no horizons at all. Simple calculations indicate that for q < q c , where q c is a critical 
value of q, it has both the inner and outer horizon, whereas for q > q c the metric potential gtt{r) has no real roots 
for r > 0. For q = q c the event and inner horizons merge and the configuration becomes extremal. It should be noted 
that the extremal ABGB geometry is distinct from, say the geometry of the extremal Reissner-Nordstrom solution. 
Indeed, because of the regularity of the former, the configurations with q > q c are not forbidden by a cosmic censor. 



A. Nonextreme black holes 



In order to determine the location of the inner and outer horizon we shall analyse equation g tt (r) — and construct 
the iterative solution restricting ourselves to the terms linear in e. As the solution can be written as 

r± = Tq + sr± , (66) 
one has to solve the simple system of two equations. First of them 

2M (rp>) 

1 - (67) 

admits solutions expressible in terms of the Lambert special functions. Indeed, it has been demonstrated in Ref. j33| 
that — Tq /M are given by 
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where W+ (s) = W(Q, s) and W- (s) = W{— 1, s) are the only real branches of the Lambert functions |48|. On the 
other hand, r± satisfies algebraic equation which can be easily solved. After some manipulations one has 

A ±] = V7TTV (69) 



1 - 2M 



Now, inserting the zeroth-order solutions (|68|l into Eq. I|69|) one can easily determine rj[ , and the result can be 
further simplified with the substitution 

- ^y, (70) 

where £ (defined as in Eq. 142(1 ') is calculated at x = xq. Specifically, denominator in the right hand side of Eq. I|69l) 
reduces to the simple expression: 

21// ' 1 : ' ' 



l^-sg- Pi) 



4£ 

To avoid unnecessary proliferation of long formulas we shall not present the results for r± here as they could be 

readily obtained by a direct substitution of the equation l|68|) and the expression describing M\ at into Eq. H69[l . 

Although the complexity of the expression describing r + makes its direct analytical application rather intricate, 
one can obtain interesting and important information analyzing limiting cases. Below we derive expansion valid for 
small q. The extremality limit will be analyzed in the next subsection. Expanding r + and collecting the terms with 
like powers of q, after massive simplifications, one has 



'■- - il, -T^-T^-^+-+ra*+(ra* + B?)^ 



< 71 . " Af + (™ k + « p. 1 ),« + ... (72) 



\320M 160M^J H \6A0M 640AT 6AM 
Analogous expression for the event horizon of the Reissner-Nordstrom solution in quadratic gravity reads 

r+ = 2M - )-Mq 2 - \Mq A - ^-Mq 6 + ... 
2 8 16 

+ (^' 3 + W 4 + ^ e + ->- < 73 » 

Inspection of (|72ll and l|73l) shows that for = (or, equivalently, k = 0) both expansion coincide up to q A terms. 

One of the most important characteristics of the black hole is its temperature. To investigate how it is modified by 
the quadratic terms, we employ a general expression (|60|l . which, in the present context, can be rewritten in the form 



1 



All r 



1 2e 



T H = — lim e-**! Si + 8nrTi . 



(74) 



Further, expanding Tu , in the auxiliary parameter, collecting the terms with like powers of e and linearizing the thus 
obtained result one has 



T H = T + sAT, 

where Tq coincides with the temperature of the nonextremal ABGB black hole 



and AT is given by 



To 



ar I 



1 

An 



1 



Q___ 

r Q Mrl 



am) 



1 



2nr^ r 



27rro 



S t t +T iP 1 (r ) + 2r T t t 



(75) 



(76) 



(77) 
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Here tq = r^ + \ whereas Tfff\ t and II^ are given by the first and the second term in the right hand side of Eq. I|32|l . 
respectively. 

Of course, the general expression for temperature is also too complicated to be displayed here. We shall, therefore, 
calculate its limiting behaviour for small q precisely in the same manner as it has been done with the location of the 
event horizon r+ . As the temperature of the extremal configuration is zero, this condition can be used in establishing 
relation between Q and M. This will be done in the next subsection. The expansion of the temperature of the 
nonextremal black hole reads 

J-H = 



8ttM 128itM 768ttM 

k q 2 a . q A /33, 



whereas for the Reissne-Norstrom geometry in quadratic gravity one has 

rp \ 9 4 9 6 | | ( g 4 , 3g 6 f 

H 8nM 128nM 128fM "' ^320M 512M '"' ' 



B. Extreme black holes 



In this section we shall investigate the important issue of extremal black holes. When r + and r_ of the classical 
ABGB black hole merge into the one degenerate horizon one has the extremal configuration characterized by vanishing 
surface gravity (Hawking temperature) . The geometry of the vicinity of the degenerate horizon belongs to the AdS2 x S 2 
class with different modules of curvatures of subspaces jl^. One of the peculiarities of the general extremal solution 
is infinite proper distance between two points, one of which resides on the event horizon. The degeneracy means 
that in the vicinity of the horizon the leading behaviour of the metric potential gu is of the form g u ~ (r — r + ) 2 , or 
equivalently 

g tt (r+) = g' a (r+) = 0. (80) 
Combining these equations one obtains simple algebraic equation that could be easily solved 

*+ = ( 81 ) 

It can be demonstrated that the parameters q and x x of the classical extremal ABGB black hole can be expressed in 
terms of the principal branch of the Lambert function evaluated at e _1 . Employing definition of W+ one has 



q c = 2^W+ (e- 1 ) w 1.0554, (82) 
whereas the location of the degenerate horizon is given by 



1 + W+ (e- 1 ) 



0.8712. (83) 



Here and below x (with or without sub or superscripts) denotes appropriate radial coordinate divided by M.. 

Now, we shall examine modifications caused by quadratic gravity. Lengthy but routine calculations show that the 
extremal canfiguration is possible for 

q 2 = q 2 + 6, (84) 

where the correction S is given by 



^(^ + ^)Li,H + M t2 (85 ) 

j=0 



with 



h » ] - - i ^o5 ! ( 516 - mw + lW - 5w3 - w ^ - ^ 2 ~ 252 °) ' (86) 
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h{ a) = - (1 ± W J (119 -Aw- 6w 2 - 4w 3 - w A ) , (87) 



h { 2 a) = -— (1 + to) (23 -Aw- 6w 2 - 4to 3 - w 4 ) , 

— = —^- 2 h ( t = jt 1 -^!^ = 4 Q) - (89) 
{1 + wf 3 {1 + wf (l + w) 5 6 w 3 

^ - JJ ^- ( 356 + 21w + 21w2 - 5 " 3 - + ^bl6€? ( 1680 + 98 " 2 - 31 " 4 ) ' (90) 

= J]t W l (257 + 28to + 2w 2 - 12w 3 - 3w 4 ) , (91) 
1 240w 3 v ' 



hi P) = — — — ~ (1 + to) (19 + 4to + 2w 2 - Aw 3 - w 4 ) , (92) 
low 5 



hf ] = -^(l + wf (l-2w-w 2 ) , hf ] = (1 + w) (1 + 6w + 3w 2 ) , (93) 



1 + tu 5 6 2w 



and 



W+ (e- 1 ) . (95) 



The term containing b% (= — /i), i. e., the last term in right hand side of Eq. I|85|l has been singled out for convenience. 
For the extreme black hole fi reads 

/« = „. , cr, (96) 
64A4-U; 3 

where cr is defined through the relation 141(1 . The location of the degenerate horizon is given by 

x = x c + A, (97) 

where A can be compactly expressed in the form 

Sw 2a + 3 , „.. » 4to(to — 1) , ,„„. 

A = 2 + ttt— w + 3 (to- 1 - — \ -U 2 , 98 

(1 + to) 2 1&wM 2X A ; Af(l + w) 2 v ; 

with <5 given by Eq. 1)85(1 . As both 5 and A depend on the particular value of the Lambert function at e _1 one can 
easily calculate their numerical values. Indeed, simple calculations yield 

0.81451 3.67845 
= Af 2 Al 2 " (99) 



and 



. 1.40646 . 3.88206 

A = ^ + ^ Q ' (1 ° 0) 
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for 62 = — fi, whereas with 62 = one has 



and 



0.57553 , 0.05121 



0.43288^ 1.05314 
A = -n^ (3 ~n^ a - (102) 

Note that for particular choices of the coupling parameters it is possible to reduce either x+ or q c to its general 
relativistic values. Obtained corrections can be contrasted with the analogous results calculated for the extreme 
Reissner-Nordstrom black hole in quadratic gravity: 

,1 = 1 + ^ (103) 



and 



(1 ° 4) 



V. CONCLUSION AND SUMMARY 



In this paper we have constructed perturbative solutions describing spherically symmetric and static black holes to 
the coupled equations of fourth-order gravity with the source term given by the stress- energy tensor of the nonlinear 
electrodynamics. The Lagrangian of the nonlinear field is a natural generalization of the Bronnikov's Lagrangian. 
Because of technical complexities we have restricted ourselves to the first-order corrections, i. e. the terms proportional 
to e. The obtained line element is parametrized by two integration constants and free parameters. Integration constants 
are related to the (magnetic) charge and a total mass of the system as seen by a distant observer whereas the free 
parameters are adjusted to make the resulting solution regular at the center. The regularity heavily relies on the 
form of the zeroth-order solution which coincides with the solution of the ABGB-type and the special properties of 
the polylogarithms Li^. The metric potentials thus computed enabled construction of the basic characteristics of the 
geometry and its asymptotics: location of the inner and outer horizon, Hawking temperature and the relation between 
Q and M for the extremal configurations. 

If one is interested in the solution possessing Schwarzschild asymptotics as Q — ► or intend to study the external 
region without any relations to the issue of regularity it suffices to put 62 = throughout the paper. We intentionally 
presented basic formulas in a form allowing for a different choices of the parameter 62. It should be stressed that the 
only choice leading to regularity of the solution at the center is b 2 = —fi. 

The calculations and results presented in this paper suggest some generalizations and obvious extensions. First, one 
may attempt to go beyond first order calculations in order to establish full regularity of the solution. Of course it would 
be interesting and desirable, with all conceptual limitations of the method, to demonstrate it in a nonpcrturbative 
way. However, the technical complexities may be a real obstacle in this regard. Further, it is possible to analyse 
behaviour of the test quantum fields in the ABGB background with the special emphasis put on the back reaction on 
the metric. Indeed, since the general expression describing the stress-energy tensor of quantized massive scalar, spinor 
and vector fields in the large mass limit is known |32t [50] it is possible to investigate the resultant quantum-corrected 
geometry even in the vicinity of the center of the black hole. Recently such a p rogramme has been carried out in the 
case of the quantum-corrected Reissner-Nordstrom black holes (see Refs. |5l|, |53 and the references cited therein). 
Finally, it should be noted that another important choice of the boundary conditions 

f r± if i = n 

"*•+>- {0* in>i' (105) 

which is, in turn, related to the horizon defined mass of the black hole is not considered in this paper. Here we 
remark only that parametrizing solution by the exact location of the event horizon and the charge is quite natural 
and certainly deserves further study. A lesson that follows from investigations of the corrected Reissner-Nordstrom 
solution is that some relations, as for example these for the extremal configurations, are simpler and physically more 
transparent. We intend to return to these problems elsewhere. 
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APPENDIX 



In this appendix we sketch the method of evaluating the indefinite integrals appearing in the solution of the equation 
(13811 . After substitution of auxiliary variable u — x^ 1 we have to find an expression for the integral of the form: 

J u p (ta,nhu) s du, (A.l) 

with natural p, s. We show that when s > 1 the above integral can be reduced to the integral with s = 1 and then we 
will express the latter by a sum containing polylogarithm functions. At the beginning we use the formulas (1.4.22.3) 
and (1.4.22.4) from the very extensive tables of integrals [53|: 



u p (tnnhuy n du 



^ + |(-ir'C)/ 



(coshu) 



2n-2k dU ' 



(A.2) 



u p {tanhuy n+L du 



sP tanh u du + 



(-1) 



n+k 



k=0 



2n-2k \k 



(coshu) 



2n-2fc 



-P 



(cosh u) 



In— 2k 



du 



(A.3) 



Evaluation of integrals containing even powers 2n — 2k > 2 of sech function can be further reduced to the case with 
the square of sech in the integrand by repeated use of the formula (1.4.24.1) of [53^: 



-du 



pu 



p-i 



(coshu)" 3 (q — l)(q — 2)(coshu) <? 2 (q — 1) sinhu(coshu) 9 1 



P(P - 1) 



du 



q-2 



(q-l)(q-2)J (coshu)?- 2 q-lj (coshw)?- 2 



rdu. 



As a simple integration by parts yields: 



-du = u p 



tanh?i — p J u p 1 tanhudu, 



(A.4) 



(A.5) 



(cosh u) 2 

we clearly see that to accomplish our calculation we have to find the integral: 

i p t&nhudu (A. 6) 

which requires special treatment. It will be proven below that it is expressible in terms of the polylogarithm functions 
Li„(z),n = 1,2, ...,p + 1, with a properly chosen argument z. There is considerable literature on analytical and 
numerical properties of these functions j45L being a generalization of Euler's dilogarithm. In the unit circle 
polylogarithm of integral order m > 1 can be defined by 



k=l 



From this series we easily derive the important recurrence relation: 



dz z 



which will also be used in the integral form: 



Li n (z) 



Lin-i(z) 



dz + C. 



Using this property and the integral representation of dilogarithm 

JO 1 



(A.7) 



(A. 



(A.9) 



(A.10) 
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one can extend the definition of polylogarithms to low orders: 

Lii(z) = - ln(l - z), Li (z) 



1 - z 



(A.ll) 



Starting from the recurrence relation in the integral form we repeatedly integrate by parts using the expression for 
the derivative of a polylogarithm (IA.8|) at every step: 

Li„(z) = C+ (In z)' U n -x(z)dz = ... = 



C + lnzLi„_i(z) - i(lnz) 2 Li n _ 2 (20 + ^(ln z) 3 Li„_ 3 (z) + ... 

(_1\n-l /I \n— 1 (-'\) n - 1 

1 ' -(lnz)"- 2 Li 2 (z)+ ) ' — (lnz)"- 1 ln(l - z) + 1 ' 



(n-2)! 



(n-1)! 



(n-1)! 



(lnz) r 



1 - z 



-dz, 



(A.12) 



where in the final step we have taken into account the form of Lii(z). If we now set z — —e 2u we recover the required 
integral in the last term: 



(-1) T 



(lnz) r 



1 - z 



-dz 



-du 



— I u n 1 tanh u du 

n 



(A.13) 



After rearrangement of the sum in the identity l|A.12fl and change n— 1 — > p we get the expression 

/ p+i p | 

u p t<mhudu = C + — - + u p ln(l + e~ 2 ") - V . , V ' — u p - k U k+1 (-e~ 21 
p + i f-^2 k ( P -ky. 



(A.14) 



The hyperbolic functions and monomials in u variable arising in our calculation can also be cast in the unified way 
when the special forms of low-order polylogarithms IjA.lljl are used. 

The case studied here falls into a wide range of expressions which can be integrated with the help of polylogarithms. 
Some occurrences of these functions in physical problems are mentioned in the Levin's book |45j . It is remarkable 
that they arise in Feynman diagrams integrals, in particular in computation of quantum electrodynamics corrections 
to the electron gyromagnetic ratio |55j . 
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